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ABSTRACT 

Two  independent  methods  are  described  for  calculating  the 
multiple  scattering  distribution  for  projected  angle  scattering 
resulting  when  very  high  energy  charged  particles  traverse  a 
thick  scatterer.  The  single  scattering  law  for  projected  angle 
scattering  is  taken  to  be  the  Rutherford  scattering  law  for  pro- 
jected angle  scattering  modified  at  small  angles  by  electron 
shielding  and  at  large  angles  by  a nuclear  form  factor 
which  gives  the  effect  of  the  finite  nuclear  size.  The  calcu- 
lations can  be  carried  through  for  any  reasonable  choice  of  X 
and  have  been  carried  through  for  two  suggested  choices  of 
for  the  examples  of  fast  ti-moson  scattering  in  2 cm  and  5 cm 
thick  lead  slabs,  with  good  agreement  for  the  two  methods  of 
calculation.  The  results  are  compared  with  the  theories  of 
Moliere  and  Olbert, 
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SECTION  1.  INTRODUCTION 

Several  multiple  scattering  theories  have  been  published^ 

^ E,  J,  Williams,  Proc.  Roy,  Soc,  AI69.  531  (1939); 

S,  Ooudsmit  and  J.  L,  Saunderson,  Phys,  Rev,  £8,  36  (1940); 

H,  S-,  Snyder  and  W,  T,  Scott,  Phys,  Rev,  76.  220  (1949); 

H,  W,  Lewis,  Phys,  Rev,  78,  526  (1950), 

^ G,  Moliere,  Z,  Naturforah  2a,  133  (1947);  1%,  ?8  (1948) 

^ S,  Olbert,  Phys,  Rev,  319  (1952); 

M,  Annis,  H,  S,  Bridge,  and  S,  Olbert,  Phys,  Rev,  89.  1216  (1953) 

H,  A,  Bethe,  Phys,  Rev,  89,  1256  (1953) 

which  are  concerned  with  the  angular  distribution  of  particles 
passing,  with  no  substantial  loss  of  energy,  through  a thick 
material  for  which  the  single  scattering  law  is  the  Rutherford 
cross  section  modified  at  small  angles  due  to  electron  screening. 

In  this  paper  we  shall  present  two  methods  by  which  a 
multiple  scattering  distribution  for  projected  angle  scattering 
can  be  obtained  from  more  general  single  scattering  laws  than 
those  already  treated.  We  shall  be  concerned  in  particular  with 
single  scattering  which  is  Rutherford  for  a wide  middle  region, 
but  is  modified  both  at  large  as  well  as  at  small  angles.  In 
scattering  from  atoms  the  modification  at  small  angles  is  due 
to  electron  screening,  while  at  large  angles  there  can  be  modi- 
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fications  due  to  additional  forces  { as  in  the  scattering  of 
'll  mesons  from  carbon  ) or  due  to  the  finite  extension  of  the 
nucleus  and  the  subsequent  deviations  from  Coulomb’s  law  inside 
the  nucleus.  We  consider  only  the  case  where  relatively  small 
are  important  so  the  approximation  sin  'K,  tan  is 

appropriate , 

There  has  been  a good  deal  of  experimental  interest  recently^'^ 

^ E,  Amalc’i  and  G,  PidecarOj,  Nuovo  Cimento,  535  (1950); 

W,  L,  Whittemore  and  R,  P,  Shutt,  Phys,  Rev,  88,  I3I2  (1952); 

E,  P,  George,  J,  L,  Redding,  and  P,  T,  Trent,  Proc,  Phys,  Soc, 

a66,  533  (1953); 

B,  Leontic  and  A,  W,  Wolfendale,  Phil,  J'ag,  1091  (1953)* 

^ M,  L,  T,  Kanangara  and  G,  S,  ^hrlkantla,  Phil,  Mag,  44»1091  (1953)* 

in  the  multiple  scattering  of  jx-mesons.  Experimental  results,  after 

\ 2 1 

being  compared  with  the  results  of  the  Mollere  and  Olbort-^ 
multiple  scattering  theories  have  been  Interpreted  as  indicating 
the  existence  of  an  anomolous  p,-rae son-nuclear  Interaction,  Among 
the  many  difficulties  arising  in  the  interpretation  of  these  experi- 
ments (see  Appendix  B),  one  of  the  most  obvious  seemed  to  be  the 
absence  of  a reliable  estimate  of  the  expected  multiple  scattering 
distribution  from  extended  nuclei. 

In  the  Mollere  multiple  scattering  theory  the  nucleus  is 
treated  as  a point  charge.  The  single  scattering  cross  section 
is  taken  to  be  the  Rutherford  cross  section  modified,  at  small 
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angles  only,  due  to  electron  screening.  This  gives  more  scat- 
tering for  large  angles  than  would  be  expected  from  extended 
nuclei.  In  the  Olbert  theory  an  attempt  is  made  to  estimate 
the  effect  of  the  nuclear  extension  by  multiplying  the  Moliere 
projected  angle  single  scattering  law  by  a step  function  which 
cuts  off  all  scattering  beyond  a given  projected  angle  ~ 

This,  however,  gives  a very  great  underestimate  of  the  multiple 
scattering  for  angles  larger  than  CJ?,  , as  for  large  angles 
Olbert 's  distribution  falls  off  as  a Gaussian  which  soon  is 
much  smaller  than  even  the  coherent  part  of  the  single  scat- 
tering law. 

In  what  follows  we  shall  outline  two  distinct  procedures 
by  which  a multiple  scattering  distribution  can  be  obtained 
from  given  single  scattering  laws.  We  deal  with  projected 
angle  scattering  as  this  is  the  usual  experimental  pareuneter. 

Our  second  method  could  be  extended  to  include  the  total  scat- 
tering angle.  Sections  2 and  3 describe  the  two  methods  we 
have  developed  for  dealing  with  multiple  scattering  problems. 

The  results  of  these  two  sections  have  been  used  to  calculate 
the  expected  multiple  scattering  distribution  of  relativistic 
(i-mesons  (cp  = 1 Bev)  passing  through  2 and  5 cm  of  lead.  Agree- 
ment between  the  two  methods  is  excellent.  Appendix  A contains 
a discussion  of  the  single  scattering  cross  section  used  for  the 
above  calculation.  Appendix  B gives  a review  of  the  experimental 
situation  with  regard  to  p-meson  scattering. 
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SECTION  2.  DESCRIPTION  OF  THE  FIRST  METHOD  OF  COMPUTING  MC^) 

The  first  method  to  be  described  consists  of  folding  together 
several  simpler  distributions  in  a manner  somewhat  analogous  to 
the  actual  effect  of  successive  layers  of  the  scatterer.  Although 
this  method  is  simple  in  principle,  it  seemed  on  first  Inspection 
that  it  would  be  extremely  tedious  to  carry  through  such  folding 
operations.  This  did  not  prove  to  be  the  case  when  the  computing 
techniques  described  below  vrere  used,  and  we  give  a rather  detailed 
description  of  the  procedure  for  this  reason. 

The  method  is  based  on  the  observation  that  if  , 

the  probability  of  a single  scattering  through  to  to 

one  side  in  projected  angle,  is  given,  then  the  multiple  scat- 
tering distribution  M(<p  ) is  completely  defined,  {Note  that  f[C^) 
is  not  normalized  to  \mi ty,  but  integration  over  all  angles  gives 
the  average  number  of  single  scatterings  in  traversing  the  sample.) 
Furthermore , if  the  actual  scattering  slab  were  replaced  by  a 
series  of  consecutive  slabs  A,  B,  C having  single  scattering 
laws  fg(cp), , where 

f(  (p  ) = f^(  f ) + 

then  the  same  m.ultiplo  scattering  distribution  results  on 

traversing  all  of  the  slabs  in  series.  If  

are  the  separate  multiple  scattering  distributions  for  A,  B, , 

then  M(  Cp)  results  on  folding  Mg, together. 
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For  the  case  of  two  components 


we  require  that  only  small  <f  are  Important  so  sin  f » tan  y f , 
and  effectively  consider  <f  as  ranRlng  from  - 00  to  + • 

The  distribution  law  for  small  angles  la  Just  the  Rutherford 
soatterlng  law  modified  due  to  electron  shielding.  Following 
Moliere  we  reproaent  this  as 

f’ 

where  f Is  the  projected  angle  and  Is  the  screening  angle 

CP„  = [1.13  1 3.76  lZ/137/  )= 

and  Q = M (N  t/A)  (Z  e=/pv)^  Here  (H  t/A)  gives  the  scattorer 
thickness  in  atoms/cm^,  p Is  the  momentum  of  the  Incoming 
particle,  % = Vpi  v Is  the  Incoming  velocity:  a is  the  "Permi- 

Thomas  radius"  of  the  atom  = 1.67  x ^ 1= 

atomic  number  of  the  scattering  material  (the  incident  particle 
is  taken  to  be  singly  charged)  and  m^  Is  the  electron  mass. 

The  modification  In  the  above  distribution  law  at  larger 
angles  Is  given  by  multiplying  f'(?)  above  by  %i<f  / f.  'I  »<> 
give 

*^11  (iiscussed  in  Appendix  A, 
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The  trick  of  this  method  consists  of  selecting  some  angle 
which  Is  a little  smaller  than  the  R,  M,  S,  width  of  the 
gausslaji  which  approximates  M(  ^ ) at  small  angles,  ) is 

then  set  equal  to  f ( ) for  | 0 \ f\  ^ t 

while  fg(  CJ  ) = 0 for  ( <f  ^ and  = f ( <p  ) for  I I ^ • The 

multiple  scattering  due  to  f^  Is  given  with  good  accuracy  by  the 
gaussian 

M^(  f ) - (2l.«p=>  (S) 


where 


(6) 


since  1 forf^lS^  • Also  a reasonably 

thick  scatterer  so  we  can  simplify  the  resulting  expression 

,9yfM  r-  


< 9% 


^/2j  + ^In  (x+^Jir+ 

o 

Q [ in  (2  9’/fJ  - ij 


so 

(7) 


It  is  convenient  to  use  the  parameter  ^ = <p / <po  and  re- 
place f(<p),  Q,  <p^  , and  <p  by  the  equivalent  quantities  g(y); 

B,  y^,  and  j' , For  the  case  of  a n-raeson  with  cp  = 1 Bev  and  a 
2 cm  lead  scatterer,  = 0,030)4.  radians  = l,7l4-°,  B = 0,126,  and 
y^  = 9,6  X 10”^,  For  the  choice  y^  = 0,^  we  obtain  ^y^=  0,75» 
showing  that  this  is  a satisfactory  choice  for  y* , Then,  for 
this  case, 

n /o  —2  /*,  r^r\ 

(8) 


M^l(y)  = (1.50  e"y 
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whore  the  subscript  1 on  emphasizes  that  this  method  ^ives 

the  distribution  in  terms  of  y = (The  distribution  Tvinction 

MgCx)  of  the  next  section  is  given  in  terms  of  the  angle 
characteristic  of  the  Moliere  theory.) 

The  multiple  scattering  distribution  Mp]^(y)  corre spot. ding 
to  single  scattering  angles^  <f  ~ obtained  u"ing  the 

following  considerations.  If  we  chose,  not  the  full  scatterer 
thickness,  but  some  sufficiently  small  fraction  , then  the 
single  scattering  law  g(y)  for  |y)>y  w'ould  have  essentially 
unit  weighting  for  no  scatterings  at  all,  and  very  small  vreight- 
ings  for  all  | y|  > y* . The  corresponding  multiple  scattering 
distribution  would  then  be  Identical  to  the  single  scattering 
law  for  |y/  > y*  since  there  is  negligible  likelihood  of  more 
than  one  scattering.  If  this  distribution  is  folded  together 
with  Itself  the  multiple  scattering  is  obtained  for  fraction  2oc. 
of  the  total  thickness.  This  distribution  differs  from  the 
single  scattering  law  by  double  scattering  terms  proportional  to 
^ , so  the  fractional  deviation  is  proportional  to  oC<<  1,  If 
the  multiple  scattering  law  for  thickness  2<K.  is  folded  together 
with  Itself  the  multiple  scattering  law  for  thickness  results. 

For  I y|  > y*  this  distribution  has  a fractional  deviation  from 
double  that  for  ZK  of  an  amount  double  that  of  the  previous 
process  (always  neglecting  higher  order  effects).  Thus  the  net 
fractional  correction  from  the  single  scattering  law  is  3 times 
that  obtained  in  the  first  folding.  Subsequent  foldings  for  8o<: , 
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l6«i  , 32*<  , etc.  have  net  fractional  corrections  (1  + 2 + l|.)» 

(1  + 2 + 4+  8),  (I  + 2 + I4.  + 8 + 16),  etc,  times  that  of  the 
folding  starting  with  thickness  , These  sums,  7»  15»  3l» 
would  be  8,  I6,  32,..,  if  we  had  started  our  folding  process  with 
an  infinitesimal  fraction  of  , since  going  from  /2  to  ©<.  gives 
1/2,  from  oC /\\.  to  oC /2.  gives  l/l}.,  etc,  for  the  extra  series  terras 
(I/2  + l/i|.  + .,,.)  = 1,  We  thus  note  that  the  net  fractional  cor- 
rection after  any  stage  of  folding  (neglecting  higher  order  ef- 
fects) is  lust  double  the  fractional  correction  obtained  by  the 
given  step  (l,e,  comparing  the  multiple  scattering  for  y ^ 0 for 
the  given  thickness  with  double  that  for  half  thickness).  We 
make  use  of  the  above  feature  to  select  as  our  starting  thickness 
a fraction  2~^  of  the  total  for  which  the  fraction  correction  of 
any  of  the  points  in  the  first  folding  is  not  too  large.  This 
fractional  correction  is  then  doubled  to  acco\int  with  good  ac- 
curacy for  the  effect  of  not  starting  the  process  with  an  in- 
finitely thin  sample.  In  the  calculations  for  2 cm  lead  and 
cp  = 1 Bev  it  was  found  to  be  suitable  to  start  with  I/8  the 
sample  thickness,  while  I/I6  the  sample  thickness  was  suitable 
for  5 cm  lead.  The  actual  folding  opex’ations  can  be  carried 
through  with  good  accuracy  by  replacing  the  continuous  fg(  (j> ) 
by  a discontinuous  function  having  valxies  only  at  regular  grid 
points.  For  the  2 cm  lead  case  where  y = 0,5,  the  region  0,5 
to  0,7  was  represented  by  0,2  g(0,6),  etc,  so  values  were  de- 
fined only  at  y^  = 0,  + 0,6,  + 0.8,  + 1,0,  with  a maximum 
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= 6,0  used  in  the  calculations.  For  the  $ cm  lead  case  we 
used  y = 1,2  and  0,1|.  Intervals  to  y = 10,0,  It  was  found  con- 
venient and  permissible  in  these  cases  to  neglect  contributions 
to  values  of  y 0 where  | yj  < y*  • (In  cases  where  these  contri- 
butions are  not  negligible  it  requires  only  slight  additional 
computing  time  to  include  them,)  The  remaining  computing  tricks 
can  best  be  Illustrated  in  terms  of  the  example  of  the  calculation 
for  the  2 cm  lead  case  where  we  started  with  I/8  of  the  total 
thickness.  Let  lumpsd  weights  given  to  the  grid 

points  0,  + 0,6,  + 0,8, , + 6,0  to  represent  the  single  scat- 
tering law  (0,2/8)  g (yj)  for  Jyjj>  with 


X°(0)=1-2Z  / (7j) 


A calculation  sheet  is  now  prepared  in  which  values  of  y^  = 


0,  0,2,  0,1|.,  , 6,0  are  entered  in  the  first  column.  The 

values  of  are  entered  in  the  second  column.  An  uncor- 
rected multiple  scattering  distribution  then  com- 


puted where 


Use  is  made  of  the  fact  that  terms  for  appear 

twice  and  we  can  write 

The  calculations  thus  proceed  as  follows.  In  the  third 
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column  the  values  of  are  multiplied  by  Lo), 

The  first  nvimber  appearing  in  the  column  (yj  = 0 here)  is  circled 
and  is  only  counted  once  in  the  subsequent  summing.  The  next 
coliunn  entry  starts  at  y.  = 1,2  and  contains  products  x.  (0»6) 

yCiXy^  - 0,6),  Again  the  flrht  term  at  y.  = 1,2  is  circled. 

The  next  column  starts  at  y.  = 1,6  and  contains  terms 

^^(yj  - 0,8),  Subsequent  columns  are  foinned  similarly  to 
produce  a triangular  array  (requiring  less  than  1 hour  of  slide 
rule  computing  time).  Contributions  from  negative  y^^  are  usually 
negligible  but  the  terms  for  the  first  few  negative  veilues  of  y^^ 
can  be  added  if  necessary.  The  rows  are  then  added  as  indicated 
by  Eq,  11  to  give  The  function 

compared  with  2 ^ ® differences  are  ad- 
ded to  to  obtain  the  corrected  multiple  scattering  distri- 
bution This  is  similarly  folded  with  itself  to  (^Ive 

and  a repetition  gives  /^^7j  ^ which  correspond^  to 
Mgj^(y).  Values  of  M^j^Cy)  (Eq,8)  are  computed  for  7^.®  0»  0^2, 

0,14.,  — — , 6,0e  The  final  M^^(y)  is  just  the  sum  of  gaussiaris 
centered  at  0,  + 0,6,  + 0,8,  -- — + 6,0,  , ’ 1 

iW.CjL)  - 2 

The  value  of  M^(7j)  ZIot  any  y^  is  computed  using  Eq,  12,  For 
yj  near  6,0  the  contributions  from  yj^  ^6,0  must  be  estimated. 
Fortunately  the  single  scattering  g(y)  and  ^(y^j)  are  decreasing 
positive  functions  of  y so  the  terms  In  Eq,  12  have  a maximum  for 
^ ^ and  the  remainder  can  be  estimated  by  noting  the  behavior 


page  eleven 


J 


/ 


naiIWi 


1- 


!■ 


R-83 


of  the  series  for  smaller  values  of  y^» 

The  above  observations  suggest  an  easy  method  for  estiifiating 
Mj^(y)  for  large  y In  terms  of  M^^(y)  and  Mg^(y)  without  actually 
carrying  through  the  process  Indicated  by  Eq.  12,  If  Mg^(y)  can 
be  approximated  by  a straight  line  on  a semi-log  plot  of  the  form 


ij  . 1..,^  ■ 1 .j 


then  Eq,  2 gives 


(13) 


2a" 


(lit) 


This  may  be  regarded  as  increasing  Mg^(yj)  by  a certain 
amount  on  the  semi-log  plot  to  produce  Alternatively, 

it  can  bo  written 

corresponding  to  a sidewise  displacement^  of  generate 

M^(y^  +A),  The  quantity  A is  obtained  by  plotting  both 
and  M„2^(y)  on  the  same  semi-log  plot  as  in  Fig,  1,  Let  Aj 
the  value  of  y at  which  has  the  same  (logarithjm^  c ) slope 

that  Mp^(y)  has  at  y = y^.  Using  Eq,  13  shows  that  ^ = <(y  )>■  /& 

50  A ~ A*  ^2,  Hero  a is  assTimed  to  be  a slowly  vnr2’’ing  function 
of  y,  (The  plot  by  comparison  can  be  made  rapidly,  using  two 
straight  edges,)  When  the  semi-log  plot  of  Mg^(y)  also  has 
curvature  so 


2 2 

%l^yk  x)  ^ MrT^(vj^)  e"^”/^e'^ 


(16) 


■ i 

. 1 
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then  the  equation  analogous  to  Eq,  1I4.  is 

- <y^> 

Mi(y^)XMBi(y^)  - <y^>  ) - <y2>  ) 

which,  for  bo  written  approximately  as 


Mi(yj)  ■JiiMgj^(yj)  e 


[1  * - <y">  )A^J 


Inspection  of  the  final  result  shows  that  thi  <?  corresponds 
again  to  a sidewise  displacement  of  by  A'^O.Sls!  to  generate 
in  the  examples  considered.  An  alternate  method  of  sidewise 
displacement  is  obtained  by  noting  that  M^^(x)  Mg^(y-x)  has  its 
maximum  at  A , about  which  it  resembles  a displaced  gaussian. 

This  suggests  that  M^(y^. ) is  generated  mainly  from  Mg^(yj^)  in  the 
region  y^^  y^  -A  so  M^^Cy)  should  be  generated  using  Mg^(y  - ^ ) 
multiplied  by  e ^ ^ , (The  last  factor  is  the  ratio 

Mai(0).)  These  methods  of  generating  the  approximate  curve  for 
M^Cy)  can  be  carried  out  rapidly  by  simple  displacements  on  the 
semi-log  plot  and  are  quite  instructive  in  giving  insight  into 
the  behavior  of  M^(y).  The  difference  between  Mg^(y)  and  gg(y) 
decreases  raoidly  if  the  choice  of  the  dividing  angle  y'  is  in- 
creased, This  is  compensated  largely  by  an  increase  in  the  width 
of  the  gaussian  M^-j^Cy),  and  thus  in  the  required  sidewise  dis- 
placement a'  to  generate  M^(y)  from  Mg^(y).  It  is  of  interest 
that  the  second  method  (of  the  next  section),  at  large  angles, 
just  folds  khe_sj,ngle  scattering  lav;  at  large  angles  with  a 
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gausslan  characteristic  of  the  Mollere  theory,  which  is  somewhat 
wider  than  for  the  method  <>r  selecting  y*  discussed  above. 

Pigs,  1,  2,  and  4.  show  the  characteristic  functions  when 
only  elastic  processes  are  Included,  It  is  seen  that  M^(y)  re- 
mains considerably  greater  than  gg  in  a region  where  the  Mollere 
curve  is  almost  the  same  as  the  point  nucleus  single  scattering 
curve.  This  Is  consistent  with  the  above  discussion  since  the 
logarithmic  slope  of  the  point  nucleus  g(y)  is  much  smaller  than 
the  extended  nucleus  g(y)  so  the  required  sidewise  displacement 
1s  less,  and  the  result  of  a given  sidewise  displacement  is  also 
less  (using  Eq,  I??,),  When  inelastic  scattering  is  included. 

Figs,  3 and  5»  the  curves  approach  those  for  a point  nucleus, 
decreased  by  a factor  Z"^  at  very  large  y. 

It  is  instructive  to  make  a further  comparison  of  this  method 
and  the  second  method  (of  the  next  section)  for  y ;>  4»  say.  This 
method  folds  together  the  small  angle  multiple  scattering  gausslan 
with  the  large  angle  multiple  scattering  distribution,  while  the 
second  method  essentially  folds  together  a wider  gausslan  and  the 
law  for  single  scattering.  The  second  method  uses  an  expansion 

In  terms  of  the  parameter  (4G)  ^ described  below,  neglecting 

—2 

terms  of  order  (4^^)  and  beyond.  This  very  roughly  corresponds 
to  neglecting  multiple  large  angle  scattering,  so  the  first 
method  would  be  expected  to  be  more  reliable  In  cases  where  such 
effects  are  not  negligible. 

As  a final  point  of  Interest,  we  note  that  the  curves  given 
In  terms  of  y apply  for  all  relativistic  momenta  where  1, 
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In  the  example  considered,  7 Is  in  units  of  1.714-  Bev  degrees* 

Thus  multiplying  the  absissa,  and  dividing  the  ordinate  by  1. 74- 
gives  the  distribution  per  Bev  degree  unit  of  ( c pf)*  This  is 
shown  by  noting  that  f ( 'f  )d  f a g(y)dy  so 

g(y)  = I (y^  + where 

o ” [1.13  + 3.76  (Z/I37(i>  )^~\  and 

^oe 

p p O 

r„^  = e /me  fi  2*82  x 10  -^cm  are  independent  of  p for 
000  ^ 

1. 

SECTION  3*  MODIFIED  MOLIERE  THEORY 
A*  Methods  and  Notation 

To  introduce  the  mathematical  methods  suad  notation  used  in 
this  -section  we  review  the  derivation  of  the  general  expression 
for  the  projected  angle  multiple  scattering  distribution  for  an 
arbitrary  single  scattering  cross  section.  The  derivation  fol- 
lows  that  of  Moliere  and  Clbert-^  and  applies  when  only  relatively 
small  angles  are  important. 

If  f ( )d  f is  the  probability  that  an  incident  particle 
undergo  a single  scattering  through  the  projected  angle  f to  f + d f 
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In  Its  passage  through  the  scatterer,  then 

OQ 

where  a.  , Is  the  probability  that  the  particle  in 

^ -«o 

passing  through  the  material  has  exactly  n single  scatterings 
through  the  projected  angles  9r  angular  inter- 
vals d d <p2»  “ — » ^ 9n*  Is  the  probability  that  there 

be  no  other  scatterings  besides  <^2  Since  we  are 

Interested  only  in  the  final  angle  cp  we  Integrate  over  inter- 
mediate angles  and  put 


* — oO 


-I 

(19) 


where  1®  probability  that  the  particle  emerge  at 

a projected  angle  between  <f  and  + d cp  after  exactly  n scat- 
terings in  the  material.  The  factor  l/nl  occurs  because  in  the 
Integration  over  intermediate  angles  the  n4  permutations  of 
----  (p^  have  each  been  counted.  Since  Bq,  18  holds  for  Independent 
events  the  order  of  the  n scatterings  is  irrelevant. 

Now  introduce  the  Fourier  representation  of  the  single  scat- 
tering cross  section.  Let 


oO 


ic<^)  = ^ ? (20) 


-oo 


Jf 


so 


■eo 


(21) 


page  sixteen 


I’ 


“tv 


r‘.  *■' 


R-83 


(where  g( ^ ) Is  not  related  to  g(y)  of  the  previous  section)*  Put- 
ting Eq.  20  into  Eq*  gives 

d9 -ov,  Q/s  if-ff-T. 'R-.1 

P iT)<i<P  = ^e 


~oC 


(22) 


The  Integrations  over  give  delta  fimctlons  which, 

after  the  ^ ^n-1  yleld|» 

^ n J 


(23) 


The  probability  of  a final  projected  angle  <P  after  any  number 

(w)  ' 

of  scatterings  is  just  the  sum  of  the  P C^)‘^/bver  all  values 


of  n. 


p(^)Jcp  = ^E(q>)c/cp^f^e  ^ 


yixo 


(24) 


where  we  use  the  result  from  Eq*  21 

<30 

[ 


^t  = 


f(«f  ) df  = g(0) 


From  this  it  can  bo  seen  that  the  solution  of  the  multiple 
scattering  problem  is  equivalent  to  the  evstluatlon  of  the  integral 

V p 

in  Eq.  24«  Molioro  has  given  an  evaluation  of  this  integral 
in  the  case  of  a scroonod  Coulomb  field.  The  main  trick  in 
evaluating  the  integral  lies  in  the  observation  that  the  existence 
of  multiple  scattering  will  smear  out  fine  grained  irregularities 
in  the  final  distribution  so  the  high  Fourier  components  give  a 
negligible  contribution, 

Moliere  used  the  single  scattering  law  of  Eq,  3 appropriate 
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for  a point  nucleus  with  electron  shielding*  We  are  interested 
in  the  case  where  this  is  modified  by  the  nuclear  "form  factor" 
as  in  Eq,  4*  Prom  the  discussion  of  Appendix  A we 


note  that 


Cl  for  small 
~ |z~'-for  large  | ^ I (our  choice) 


with  the  rapid  change  occurlng  when  ~ 2 

Olbert’s  method  consisted  of  using  a step  function  f or  ^ 


I for  1 I ^ 

o for  I qp  1 > <Po 


(Olbert) 


As  mentioned  in  the  introduction,  this  gives  a very  large  under- 
estimate of  the  multiple  scattering  for  ^ “fo  where  the  Olbert 
function  falls  off  as  a gaussian  while,  for  any  reasonable  form 
factor,  the  multiple  scattering  distribution  lies  above  the 
single  scattering  distribution  f ( ) for  large  y,  and  thus  far 

above  the  Olbert  distribution. 


B,  First  Derivation  of  M2  (x) 


We  first  treat  the  multiple  scattering  from  extended  nuclei 


as  a correction  to  the  Moliere  theory  by  setting  ^ = [1  +•  oj , 

where  ( - l)  gives  the  correction  term.  This  is  the  method 

used  by  Olbert  for  his  step  function.  Such  a treatment  yields 
satisfactory  results  for  small  angles,  but  is  Inconvenient  for 
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large  angles  where  a second  method  (described  following  this)  Is 
preferred,  Eq,  I4.  Is  written 

f(<P)  = I- c CX-O]  . 


The  first  Integral  has  been  treated  by  Moll^re*^  and  to  a 


See  H,  A,  Bethe,  Reference  4»  ® discussion  of  this  Integral, 


sufficient  degree  of  accuracy,  la 

a ~ G!  r f ^ 


Wo  now  Introduce  the  various  parameters  typical  of  the  Mollere 
theory,  following  the  notation  of  Olbert, 

^ j 7=  (30, 

1 

and  ^ 6 +-  L/.«3y3^-f-  3.76  (Z//37)’^  J 


la  Independent  of  the  momontvim  in  the  relativistic  region. 
Putting  Eq,  29  In  28,  and  using  3O  gives 

Where  r . f (e  '^  l)  0<  ~ 

-h 
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Inserting  this  result  into  Eq*  214-  we  obtain  in  terms  of  the 
parameter  ^ ^ ‘ T ‘V"'  ^ ll 


We  use  the  fact  that  G is  a reasonably  large  number  (for  moderate 
foil  thickness  G varies  from  about  7 to  1$)  and  expand  the  portion 
of  the  exponential  which  is  nraltiplied  by  (2G)  • This  gives, 

to  first  order  in  (2G) 


After  performing  the  integrations  over  this  become-s 

'Trr  ^ 


“horo  ^ TT/^'/xll  + -(X-x'J 

i- e - J 

Kw  = vf! 


=-|rP^e 


f’(X;«0)  is  the  Moli^re  function.  The  integral  K(x)  which  oc- 
curs in  Eq,  314.  must  be  evaluated  n\imerically  for  a given  3^  /Xo ) , 
Eq.  3I4.  is  seen  to  have  the  form  that  was  desired  originally.  The 
integral  which  contains  the  effect  of  the  nuclear  extension  oc- 
curs as  a correction  to  which  is  the 

Moliere  multiple  scattering  distribution  for-  a point  nucleus  in 
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terms  of  the  projected  angle  parameter  x.  This  corresponds  to 

•»  1 for  all  X,  If  ^(x/x^)  la  the  Olbert  step  function,  Eq,26, 

the  correction  term  becomes  the  Olbert  correction  function  K(xjx^) 

if  use  is  made  of  the  fact  that  x^  > > x^* 

o m 

Using  Eq,  34-  the  multiple  scattering  distribution  can  be 
detenulned  for  a particular  However,  Eq,  34  becomes 

inconvenient  for  large  values  of  x,  because  in  that  case  K(x) 
becomes  almost  equal  to  f ^ (x; oo)  and  the  difference  between  two 
large  numbers  must  be  used  to  give  a small  one.  For  large  values 
of  X therefore  it  la  better  to  treat  the  modified  cross  section 
directly,  rather  than  as  a correction  to  the  Moliere  distribution. 
This  is  done  as  below. 


C,  SECOND  DERIVATION  OP  M^CX), 

Consider  again  the  Fourier  transform  of  the  single  scattering 


law. 


cO 


3(0- 

— 


(35) 


Introducing  the  MolleVe  parameters  x,  and  0 as  defined 
Eq,  30  we  get 

■>?  \ ^ J-  ^ - 


and 


- ^ f' 


^ i'>X^  yO-C'n) 


( 36a) 


page  twenty-one 


R-83 


Eq*  36  of  course  la  the  Moliere-Pourier  transform  if  {x* /x^)  = !• 
In  the  case  of  the  point  nucleus  this  Integral  is  evaluated  through 
the  observation  that  there  exists  an  angle  x'  = R at  which  the 
integral  can  be  split  such  that  but  ^ where 

la  the  frequency  in  the  neighborhood  of  the  main  Fourier  compon- 
ents,^ Then  in  the  integration  up  to  H.  the  frctor  (cos  '*]  x’  - 1) 
can  be  put  equal  to  -{"v^  x*  ) /2  and,  in  the  integration  from 

R to  , (x*^  + becomes  x'*"^.  Both  integrations  can 

be  performed,  and  the  splitting  point  R cancels  up  to  terms  of 
2 

the  order  R $ which  are  very  small. 

In  the  case  with  which  we  deal  the  argument  is  modified  as 
follows.  Again  we  split  the  integral  at  an  angle  R where  x « 

X*  SI  R‘^'^'^^,  • Now  we  make  use  of  the  important  property  of 

the  form  factor  given  in  Eq,  25  that  % {x^/x^)x.l  for  small 
values  of  the  argument  x'/x^.  In  particular,  for  x'  = I/I4.  our 
chosen  fonri  of  (Appendix  A)  gives  {I/I4.)  ~ 0,96  for  the 

2 cm  lead  case,  and  0,93  I'or  6he  5 cm  lead  case.  The  integral 
up  to  R is  thus  the  same  as  in  the  case  of  the  point  nucleus, 

(37) 

Performing  the  first  Integration  and  using  the  fact  that  » 

Eq,  37  becomes 

S (38) 

where  ~ 
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It  is  understood  that  K must  be  much  smaller  than  1»  If 


now  ) Is  put  into  Eq«  3^0  and  the  terms  multiplied  by  (20) 

o.SIc'*' 

are  expanded  In  the  same  way  as  beforOjf^the  integrations,  the 
following  result  is  obtained  for  the  first  order  in  (20)“^, 

^ ^ N'g  i ^ ( 39 ) 

where  and 


-1 


-X 


TUX)  = e + e- 


iQ. 


In  its  present  form  Eq*  39  is  convenient  for  calculation  only 

for  large  values  of  x(x—  I4.  in  the  5 lead  case)  because  in  this 
-x^ 

case  e is  so  small  that  the  precise  value  of  n does  not  influ- 
ence the  result* 

However,  for  smaller  values  of  x the  exact  value  of  K is  of 
importance.  In  order  that  the  expansion  of  (cos  ''-jx’-l)  in  Eq,  37 
be  correct  U.  must  be  much  smaller  than  1,  However,  it  is  in- 
convenient to  evaluate  the  final  Integral  of  Eq,  39  numerically 

0 

for  very  small  values  of  K. 


° We  can  see  more  clearly  what  error  is  Introduced  if  K is  al- 
lowed to  become  large  by  doing  the  following,  Expgmd  T(x, X ) 

of  Eq,  39  in  powers  of  2xX  , Then  neglecting  terms  of  order 

2 

or  (xX  )^  gives  T(x,X  )-  2 X ^ e”^  (2x^-1) 
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Integrating  N(  X)  from  K/  to  L and  letting  ^ /x  ) = 1 

In  this  integral  we  obtain 

N(LjX)J 

But  this  is  Just  the  same  result  obtained  by  letting  K-  = L in  the 
Fourier  transform  s(  ) of  Eq,  37*  Thus  allowing  to  become 
large  introduces  the  same  error  as  would  be  Introduced  if  the 
final  result  were  expanded  in  powers  of  2xA» 


To  overcome  this  difficulty  one  csin  use  the  property  of  the 

form  factor  that  ^ for  values  of  A ™*ch  larger  than 

7?;'  . Then  / A'^TC^A)^>  can  be  evaluated  analytically 

where  L Tj^ig  makes  the  numerical  integration  J A 

L 

practicable.  In  doing  this  it  will  be  seen  that  the  dependence 
upon  vanishes  and  a convenient  and  accurate  expression  for 
Mp(  X.  ) is  “’obtained.  We  observe  that  if  L '^1/14.  than  A can  be 

if 

considered  small  even  if  2 A is  not.  Thus  we  can  expand  e’ 
in  T(  Xj  A ), 

_ V Ax  - 0 ^ ’‘(1-  A’+ 


T(x,A) 


(i;0) 


taking  only  terms  up  to  A , We  evaluate 


T 


J XJ  ' 


K- 


This,  after  some  labor,  «sA  neglecting  terms  in  ^ or  higher,  gives 
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T ^ '2.  L»  X _ l2L_  s', A -k  (o^  ^ + 

N.  V-  \_^^v  jV  ^ 


Putting  this  result  Into  Eq,  39  gives 


where ^ / <L-t  / v 


Vw 


_ S',  *^V\  "X Vi 


('"■V.v  4 „ (2:^  4 ^3-^  + ■■ 
3 ~ X.-x  H 


In  our  earlier  coimnunication  by  the  same  title,  Leon  N,  Cooper 
and  James  Rainwater,  Phys,  Rev,  9^.  110?  (195^4-) » q(L,x)  was  given 
as  2(2x^  - 1)  q (K)  = 2(2x^  - 1)^  (K/l.26)  which  was  a svif- 
ficiently  good  approximation  in  the  example  there  considered. 


D,  EVALUATION  OF  THE  INTEGRALS 
The  integrals  K.  (*)  N(L,  x)  which  occur  in  expressions 
for  the  multiple  scattering  distributions  derived  above  csui  easily 
be  evaluated  by  numerical  means,  N(L,  x)  in  particular  can  be 
evaluated  conveniently  for  values  of  L = l/lj.  or  1/2,  suid  for  such 
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values  of  L Eq,  43  is  accurate.  These  integrals  do  not  depend 
on  the  momentum  of  the  incoming  particle  but  rather  on  , which 
/i'  1 in  the  relativistic  region.  Thus  a single  computation  is 
sufficient  for  all  relativistic  momenta.  For  different  thick- 
nesses, or  atomic  number,  the  form  factor  is  changed  only  through 
the  change  in  in  Table  1 gives  the  function 

T()C^A)Afor  various  values  of  X at  grid  spacings  of  1/4  for  A • 
Table  2 gives  the  values  of  the  function  q(L, X ) for  L = 1/4  for 
relevant  values  of  X , 

K(x)  and  N(L, X ) have  been  evaluated  numerically  for  the 
2 cm  and  5 cm  lead  cases  for  grid  spacings  of  A x = 1/2  and  1/4 
using  Weddle's  rule,^^  Comparison  of  the  results  for  the  two 

H,  Margenau  and  G,  M,  Murphy,  ’’The  Mathematics  of  Physics 
and  Chemistry**,  D,  Van  Nostrand  Company,  Inc,,  New  York, 

N,  Y,,  1943 » page  46l, 

grids  shows  the  largest  numerical  errors  occur  for  the  small 
values  of  X where  the  correction  term  is  unimportant.  For 
lareer  values  of  X.  (X^  3)  the  change  in  the  numerical  results 
for  grid  changes  from  1/2  to  1/4  is  less  than  2 percent.  In  all, 
the  errors  in  the  distribution  due  to  numerical  errors  resulting 
from  the  integration,  for  a grid  of  1/4,  appear  to  be  less  than 
1 percent. 

For  very  large  values  of  X asymptotic  formulas  can  be 
developed.  For  a form  factor  which  decreases  asymptotically  as 
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/f)-3 

X the  asymptotic  expression  is 


"M 


r 1 " ( ",  t.i)  + 1 

L 1;  J 


Such  expressions  give  reasonably  accurate  results  for  x ^ 6 (5  P®r 
cent  or  better).  The  procedure  described  near  the  end  of  section 
2 is  also  readily  adapted  to  this  method  for  a quick  graphical 
correction  procedure  to  obtain  Mglx)  from  the  single  scattering 
law.  For  real  atoms  at  the  larger  angles  incoherent  scattering 
predominates  so  that  the  form  factor  becomes  Z”^,  Then 

the  multiple  scattering  distribution  is  given  by  the  asymptotic  ex- 
pansion of  the  Moliere  distribution  multiplied  by  Z"*^,  This  is 

M2(x)dx  = ^^1  + ^+  ^+ .J 

The  results  presented  in  this  section  have  included  only 
the  first  powers  in  (2G)“^,  It  is  possible  to  obtain  the  next 
order,  but  in  view  of  the  large  xincertainty  in  the  nuclear  form 
factor  it  was  not  considered  worth\v’hile  at  present  to 

consider  this  term. 

We  wish  to  thank  Professor  Robert  Berber  for  helpful  dis- 
tKe. 

cussions  of/^single  proton  form  factor.  We  also  wish  to  thank 
Miss  Hilda  Oberthal  for  her  assistance  with  the  computations  of 
the  results  of  the  modified  Moliere  theory. 
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The  single  scattering  laws  used  in  the  calculations  were  ob- 
tained in  the  following  way*  The  Rutherford  formula,  modified  at 
small  angles  due  to  electron  shielding,  is  given  in  Eq,  3 following 
Moliere*  When  nuclear  extension  is  considered,  this  must  be 
multiplied  by  a nuclear  form  factor  wMch  contains  a 

part  representing  elastic  coherent  scattering  plus  a part  rep- 
resenting inelastic  scattering  if/fo  ) + )• 

In  principle,  if  the  nuclear  wave  functions  were  known  exactly, 
one  would  prefer  to  calculate  Pj^  by  an  exact  phase  shift  analysis 
for  spin  1/2  particles,  and  calculate  P^  by  considering  in  detail 
all  of  the  possible  final  states  of  the  scattering  system.  How- 
ever, the  nuclear  charge  distribution,  and  the  nuclear  wave 

functions  are  not  known  exactly,  and,  in  fact,  measurements  of 
C 11 

Pjj  for  fast  electrons  are  providing  valuable  information  con- 
R,  W,  Pidd,  C,  L,  Hammer,  and  E,  C,  Rsdca,  Phys,  Rev*  92* 

436  (1953); 

R,  Hofstadter,  H,  R,  Fechter,  and  J,  A,  McIntyre,  Phys*  Rev*  92* 

978  (1953); 

L.  I,  Shiff,  Phys*  Rev*  Si,  988  (1953); 

R*  Hofstadter,  B*  Hahn,  A*  W,  Knudsen,  and  J*  A*  McIntyre, 

Phys,  Rev,  9?.  512  (1954); 

D,  R,  Yennie,  D*  G*  Ravenhall,  and  R.  N,  Wilson,  Phys,  Rev*  95« 

500  (1954) 
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mesonic  x-ray  studies,  fast  electron  scattering  and  other 


Val,  L,  Pitch  and  James  Rainwater,  Phys.  Rev.  92.  769  (19^3); 

L,  N,  Cooper  and  M,  Henley,  Phys,  Rev.  92,  801  (19^3) » 

John  A,  Wheeler,  Phys,  Rev,  92.  8l2  (1953); 

D,  L.  Hill  and  K.  W,  Ford,  Phys.  Rev.  9^,  161?  and  I63O  (19514-) 
12  11 

experimental  results  * which  are  sensitive  to  the  nuclear 

P,  Bitter  Eind  H,  Feshbach,  Phys,  Rev,  92,  83?  (1953); 

B,  G,  Jancovlcl,  Phys,  Rev,  95*  369  (19514-); 

D.  C,  Peaslee,  Phys,  Rev,  9£,  717  (1951;) 

charge  distribution,  all  agree  that  the  nuclear  charge  distri- 
bution 1s  more  compact  than  had  previously  been  believed  to  be 
the  case.  For  a nuclear  charge  distribution  (r ) = for 

r R = and/^(r)  = 0 for  r>R  the  radius  R has  clear 

meaning  and  a best  match  to  the  above  experiments  gives  r^<;^l,0 
to  1,2  X 10”^^cm  for  not  too  small  A,  Elementary  considerations 
of  quantum  mechanics  show  that  such  a model  cannot  be  strictly 
correct,  and  y^(r)  must  be  a continuous  function  of  r.  In  this 
case,  although  the  different  experiments  are  not  always  sensitive 
in  the  same  manner  to  the  shape  ofy^(r),  the  equivalent  uniform 

model  r is  usually  taken  to  be  that  value  which  gives  the  same 
o 

as  for  the  non-uniform  (r) , Prior  to  late  1953  the 
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Value  of  r would  have  been  considered  to  be  "well  known”  and  of 
o 

magnitude  l.l^.  to  1,5  x 10“^^cm,  Thus  all  of  the  comparisons^*^ 
between  the  experimental  and  expected  multiple  scattering  distri- 
butions used  this  larger  value  of  r^» 

Until  recently  It  has  been  customary  to  calculate  using 
the  Born  approximation  which  gives  a linear  superposition  of 
the  scattering  amplitudes  of  the  individual  protons  so  the  scat- 
tered intensity  is  of  the  form 


where  PiT  is  the  scattering  intensity  of  a single  proton,  q = 
represents  the  vector  momentum  change  in  scattering,  r^  is  the 
position  coordinate  of  the  proton  in  the  nucleus,  and  the 
evaluation  is  for  the  grovind  state  |i  ^ of  the  nucleus.  Thus 
is  just  the  Fourier  transform  of  the  nuclear  charge 
distribution.  For  a uniform  nuclear  modelj  letting  y , 

this  gives 

y " y y^  J 

which  gives  diffraction  minima  at  y = 7*7»  etc,  Schlff^^ 

In  particular  has  considered  the  form  of  Fjj  (y)  in  Born  approxi- 
mation for  various  simple  analytic  forms  for/^(r).  We  note  that 
Fjj  1s  of  the  fonii  of  the  square  of  a real  amplitude  term,  and 
thus  gives  diffraction  minima  when  the  amplitude  changes  sign. 

For  distributions  sufficiently  peaked  in  the  center  and  with  a 
gradual  ”tailing-off ” , the  amplitude  may  not  change  sign  and 
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F^-  will  be  a smoothly  deci'easinf':  function  of  y.  It  la  readily 

Q 

seen  from  Eq,  A1  that  Pj^(O)  = 1 and,  by  expanding  the  exponentials 
and  using  inversion  symmetry  thrcupii  r = 0,  that  the  leading  term 
in  the  decrease  of  Pj^  for  small  y depends  on  <^r  ^ , In  this 

connection  it  is  interesting  to  note  that  the  widely  used 
Williams^  formula  uses  V(r)  = (Z  e^/r)  (1  - for  which 

/^(r)  = (I^y^/3)  (R/r)  and  r^^  = 3R^/2» 

charge  density  for  a uniform  nucleus  of  radius  R,  Although 
Williams  implied  that  this  distribution  closely  approximates  a 
uniform  modal,  and  it  has  been  taken  by  others^  as  corresponding 
to  the  uniform  modal,  it  actually  corresponds  to  a rather  strongly 
peaked  distribution  about  r = 0 with  <^r  ^ larger  by  a factor  of 
$/2  than  for  a uniform  distribution,  and  thus  a correspondingly 

c 

more  rapid  initial  decrease  of  for  small  y»  For  larger  y it 

Q 

gives  a uniform  decrease  of  with  y with  an  asymptotic  form 
l6/y^  for  large  y.  In  the  region  of  large  y Eq,  A2,  between 
minima,  has  a steady  decrease  with  y which  can  be  estimated  by 
neglecting  the  sin  y term  and  setting  jcos  y | "X-  1 to  give  an 
asymptotic  form  9/y^,  which  is  belovf  that  for  the  V'illiams  distri- 
bution. 

The  detailed  phase  shift  calculations  of  Yennie,  Ravenhall, 
and  Wilson,  using  various  assumody*^  (r ) , show  that  the  shape  of 

n 

P^  (y)  is  energy-dependent  and  significantly  different  from  the 
Born  approximation  value  for  high  Z materials.  In  particular, 
the  scattering  amplitude  is  a complex  number  which  circles  the 
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value  zero  in  the  complex  plane  when  "changing  sign"*  Thus  the 
"diffraction  minima"  are  largely  missing^  or  are  greatly  reduced 
in  magnitude  compared  to  the  Born  approximation  results.  This 
feature  is  also  apparent  in  the  experimental  results  and  had  led 
to  an  initial  "Born  approximation  interpretation"  thaty^  (r)  must 
resemble  an  exponential  distribution.  The  interpretation  favored 
at  the  time  of  this  writing  is  thaty^(r)  can  be  fairly  constant 
for  values  of  r containing  most  of  the  charge  distribution,  with 
a gradual  dropping  off  at  the  "surface". 

In  view  of  the  above  results,  and  with  the  consideration  that 

n 

we  wish  to  choose  a form  for  P„  which  will  not  underestimate  the 

N — — ' 

expected  multiple  scattering  due  to  the  known  electromagnetic 

Interaction  of  mu-mesons  with  protons  in  examining  the  results 

of  experiments  Investigating  possible  "anamolous  scattering", 

c • 

we  have  chosen  the  following  form  for  Pj^,  For  y = 0,  1,  2,  3 

Q 

we  choose  Pj^  = 1,00,  0,02,  0,^0^  and  0,1^  to  approximate  Eq,  A2 

after  the  effect  cf  the  first  diffraction  minimum  is  "removed", 

A smooth  curve  through  these  points  is  then  joined  smoothly  to 

= 12/y^  for  y ^ I4.,  this  being  between  the  Williams  formula 

and  the  value  obtained  above  for  a uniform  distribution.  In 

principle  this  should  be  applied  to  the  cross  section  for  total 

angle  scattering  rather  than  for  projected  angle  scattering  as 

we  do  here.  The  consequences  of  this  approximation  are  discussed 

I 

following  the  discussion  of  F^, 

The  calculation  of  the  inelastic  scattering  is  quite  dif- 
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flcult  to  perform  exactly.  Amaldi,  Pidecaro,  said  Mariani^^^  have 
E,  Amaldi,  G,  Pidecaro,  and  P,  Mariani,  Nuovo  Cimento 

553  (1950) 

carried  out  detailed  calculations  of  the  inelastic  and  elastic 
scattering  expectea,  using  a particular  independent  particle 
nuclear  model.  For  incident  particles  of  very  high  kinetic  energy 
and  cp  100  Mev),  a Born  approximation  treatment  similar  to 

that  used  in  the  theory  of  x-ray  scattering  by  atoms  should  be 
reasonably  accurate  and  gives  the  result,  as  for  x-rays, that 


^ c,  f,  A.  H,  Compton  and  S,  K,  Allison,  "X-rays  in  Theory 
and  Experiment”,  D,  Van  Nostrand  and  Co,  (1935)*  Chapter  3* 

P^  ~ Z (1  “ Fjj)  times  the  form  factor  for  the  scattering  by  a 
single  proton.  This  last  factor  must  be  included  if  the  proton 
is  not  effectively  a point  charge  (due,  say,  to  meson  cloud  ef- 
fects), Amaldi  etal^^  have  particularly  emphasized  this  point 
and  have  calculated  the  expected  effect  of  the  proton  "size”'  on 
the  basis  of  a simple  model.  Experimentally,  however,  it  seems 
that  the  proton  charge  distribution  should  be  treated  as  being 
confined  to  a surprisingly  small  volume  on  the  basis  of  electron 
scattering  experiments,^^  and,  by  Inference,  from  the  Interpre- 

J.  A,  McIntyre  and  R,  Hofstader,  Bulletin  of  the  American 

Physical  Society,  Vol,  29,  paper  II  (also  see  paper  12)  A.P.S. 
meeting  at  Seattle,  Washington,  July,  1954« 
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tation  ' of  the  experiments  on  the  neutron-electron  interaction* 


Poldy,  Phya.  Rev.  87.  675  (1952)  for  a discussion 

of  this  subject  and  for  reference  to  earlier  papers. 

The  experiments  show  that  the  proton  form  factor  is  essentially 
tuilty  for  cp  200  Mev  and  at  angles  90°,  so  we  teOce  it  to  be 
essentially  unity  for  cp tf  values  considered  in  this  paper.  For 
lai'ge  eiiuugh  value  of  cp  <p  this  factor  will  eventually  become 
important  and  require  consideration. 

If  the  single  proton  form  factor  la  set  equal  to  imlty,  the 
above  expression  for  can  be  understood  by  analogy  with  x-ray 
scattering  where  the  ’’inelastic”  scattering  corresponds  to  modi- 
fied Compton  elastic  scattering,  where  the  recoil  momentum  is 
taken  up  by  a single  (moving)  electron  rather  than  by  the  atom 
as  a whole.  The  final  states  of  importance  correspond  to  recoil 
electron  momenta  centered  about  the  photon  momentum  transfer, 
modified  by  the  initial  electron  momentum  distribution,  which 
is  given  by  the  Fourier  transform  of  the  ground  state  wave  function 
of  the  atom.  For  charged  particles  of  kinetic  energy  and  cp  >> 

100  Mev,  and  for  y >>  1,  the  final  states  of  importance  should 
be  attainable  with  moderate  energy  loss,  so  the  expression  for 
the  total  scattering  intensity  in  a given  direction  can  be  written 
with  fair  accuracy  by  summing  over-all  final  states,  keeping 
for  a given  angle,  the  same  as  in  the  elastic  case. 
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z"[F>  Fj]  = ^ ' Ii>  I 


■ ^ 


vn 


(A3) 


by  closure.  For  an  independent  particle  model  Eq,  A1  and  Eq,  A3 


lead  to  the  above  quoted  result  for  Thus  we  set 

C+z^.-Fn') 


(AI4.) 


The  remaining  point  to  be  considered  is  the  error  inti-o- 
duced  when  is  applied  to  the  law  for  projected  (rather  than 
total)  angle  scattering.  For  a do/dO  law  varying  aa  (f>  , the 

projected  angle  single  scattering  law  Is  obtained  by  multiplying 
do/dQ  by  where  cCv%)=  k/2,  I4./3,  l6/l5,  ^Tt/l6,  and  32/3^ 

for  n = l4-»  5»  7»  8»  and  9 respectively  ( CCvi)  ~ ^ j c-trcu*  ‘P 

O 

In  going  from  lO  =(4.  for  Rutherford  scattering  to  /I  = (which 
is  obtained  when  the  asymptotic  form  of  Pj^  multiplies  the  Ruther- 
ford scattering,  c(n)  is  reduced  by  the  factor  of  5/8*  This  rep- 
resents an  extreme  situation  since  more  slowly  varying, 

V/e  note  that  the  above  effect  can  approximately  be  tsdcen  into 
account  by  choosing  r^'^  10  percent  larger  than  otherwise  when 
applying  to  the  law  for  projected  angle  scattering.  In  the 

examples  we  choose  r^  = 1.1  x 10”^^  cm,  corresponding  to  r^ 

1,0  X 10”^^  cm  for  da/dO,  In  view  of  the  uncertainty  concerning 
the  exact  form  for  the  true  consider  this  approxim- 

ation to  be  adequate  for  the  present. 
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ON  MU- MESON  SCATTERING 

The  interest  of  the  authors  in  the  multiple  scattering  theory 

was  mainly  stimulated  by  the  possible  consequences  of  a strong 

anomolous  mu-meson-nucleon  interaction  on  the  interpretation  of 

12 

the  mu-mesonic  x-ray  experiments.  Cosmic  ray  experiments  on 
the  scattering  of  fast  mu-mesons  by  nuclei  suggested  that  some 
anomolous  scattering  exists,  but  the  interpretation  of  the  experi- 
ments are  not  completely  unambiguous,  and  there  is  considerable 
disagreement  on  the  magnitude  and  existence  of  the  effect.  When 
these  experiments  were  analysed,  the  expected  multiple  scattering 
distribution  for  an  extended  nucleus  was  obtained  by  various  ap- 
proximations such  as  comparing  the  experimental  results  with  the 
predictions  of  the  Moliere  and  Olbert  theories,  or  by  using  the 
Williams  theory.  In  all  cases  the  old  ’’large”  nuclear  size  was 
used  which  gives  <r  > twice  that  favored  by  recent  experiments. 
Thus  the  experiments  were  always  analysed  on  a basis  that  under- 
estimated the  coulomb  multiple  scattering.  Aside  from  the  results 
oi"  the  experiments  discussed  below,  we  note  that  the  mu-mesonic 
x-ray  results  indicate  that  any  anomolous  energy  independent 
nuclear  potential  for  the  mu-meson  can  in  its  effect  at  most  be 
equivalent  to  a slight  change  in  the  choice  of  the  nuclear  radius 
when  calculating  the  coulomb  Interaction,  This  could  not  explain 
any  significany  portion  of  the  anomolous  scattering  reported  in 
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some  of  the  experiments.  Also,  experiments  on  the  nuclear  inter- 
action  by  Annis  at  al"*^  can  be  explained  without  invoking  any 


M,  Annis,  H,  C,  Wilkins,  sind  J.  D,  Miller,  Phys,  Rev,  9^4-. 

1038  (I95i4-) 


anomolous  interaction.  Such  an  ”Einomolous  interaction**  would  have 
to  be  strongly  energy-dependent  and  thus  only  appear  strongly  at 
high  energies  to  explain  the  low  energy  experiments, 

Amaldi  and  Fidecaro^  investigated  the  large  angle  scattering 
of  fast  mu-mesons  in  the  energy  bsuids  200  Mev  to  320  Mev  and  ^320 
Mev,  using  a counter  hodoscope.  They  compared  the  large  angle 
multiple  scattering  in  iron  and  lead,  emphasizing  the  iron  re- 
sults as  far  as  anomolous  scattering  is  concerned.  About  5 x 10^ 
incident  mesons  were  counted  and  results  for  iron  were  consistent 
with  no  anomolous  interaction.  In  their  series  2 run  on  iron  214.9, 

168  particles  were  incident  and  3 scattered  particles  were  ob- 

2 

served  in  each  energy  band.  When  an  extra  200  gm/cm  of  bricks 
were 

placed  above  the  apparatus  to  decrease  the  nvunber  of  protons 
etc,,  2OI4.,  914  particles  were  incident  with  one  scattered  count 
in  the  lower  energy  band  and  none  in  the  higher  energy  band.  The 
latter  numbers  are  about  the  expected  values  for  scattered  protons, 
etc,  with  that  amount  of  filtering.  They  set  an  upper  limit  for 
"anomolous  scattering"  of  --  4*5  ^ 10  cm  /nuclean  in  the  lower 
energy  band  and  ^ 2,3  x 10“^*^cm^/nuclean.  in  the  upper  energy 
band.  These  values  assxime  Isotropic  scattering  for  the  anomolous 
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part, 

Whlttemore  and  Shutt  Investigated  the  multiple  scattering  of 
negative  mu-mesons  in  5 cm  of  lead  for  particles  having  0,3  Bev  ^ 
cp  ^ 3»1  Bsv  using  two  cloud  chambers  and  a magnet  for  momentum 
analysis.  Their  experimental  points  essentially  fall  on  the  Scott- 
Snyder  (Mollere)  curve  for  p«f^  I3  Bev  degrees^  with  one  point  at 
17  Bev  degrees  a factor  of  about  two  below  the  Mollere  curve  (but 
with  a large  statistical  uncertainty).  Prom  Pig,  5 these  points 
would  also  fall  above  our  theoretical  curve.  The  principle  dif- 
ficulty in  the  interpretation  of  these  results  is  the  question  of 
the  certainty  that  only  mu-nesons  were  Involved,  It  is  of  interest 
in  this  connection  that  the  later  experiments  were  conducted  under- 
ground, or,  in  one  case,  using  1 meter  of  lead  absorber  above  the 

cloud  chamber  to  assure  greater  beam  purity.  Also  we  should  like 

19 

to  point  out  that  measurements  on  the  scattering  of  cp  200  Mev 

John  0,  Kessler  and  Leon  M,  Lederman,  Phys,  Rev,  9k, 

689  (195Ii) 

n-mesons  on  Pb  show  differential  cross  sections  for  elastic  plus 
Inelastic  scattering  which  are  always  <v0,l  barn  or  larger.  By 
contrast,  the  elastic  scattering  experiments  using  electrons  of 
comparable  momentum  give  differential  cross  sections  which  are 
^ 10”^  barns  at  120°,  a factor  of  10^  below  the  n-meson  cross 
section.  This  contrast  is  admittedly  extreme,  but  it  emphasizes 
the  importance  of  not  underestimating  the  possible  importance  of 
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small  percentages  of  beam  contamination  in  such  experiments. 

The  remaining  experiments  used  greater  absorber  thickness  at 
sea  level,  or  operated  underground  to  minimize  beam  contamination. 
They  have  the  common  feature,  however,  that  the  momentum  distri- 
bution of  the  incident  particles  was  not  measured  directly,  but 
was  ass-umed  known  from  other  sources.  Any  error  in  the  assvimed 
known  momentum  distribution  would  tend  to  affect  all  of  these 
experiments  in  a similar  fashion.  The  fact  that  the  observed 
scattering  distributions  are  not  given  directly  as  a function  of 
p makes  comparison  of  their  final  curves  with  our  calculated 
M(f)  difficult, 

George,  Redding,  and  Trent  measured  the  multiple  scattering 
of  penetrating  cosmic  ray  particles  in  2 cm  lead  plates  at  60  m,w,e, 
vinderground  using  a counter-triggered  cloud  chamber.  Three 
experimental  arrangements  were  used  for  the  triggering  counter 
telescope  involving  0,  5»  and  10  cm  lead  below  the  cloud  chamber. 

The  particles  were  all  assumed  to  be  p,-mesons  and  the  momenta  of 
the  individual  particles  were  known  only  to  be  above  the  cutoff 
values  determined  by  the  lead  absorber  thickness.  The  analysis 
was  made  by  assuming  that  the  energy  d.''  stribution  was  flat  for 
E <<  = 12  Bev,  The  experimental  distribution  N(^)  was  com- 

pared with  one  calculated  using  a weighted  average  of  gausslan 
fimctlons  (one  for  each  energy)  in  accord  with  the  above  pre- 
scription for  the  assumed  energy  distribution.  No  anomolous 
scattering  was  observed  using  10  cm  of  lead,  but  a small  amount 
of  "anomolous  scattering"  appeared  at  larger  angles  when  0 or  5 cm 


page  thirty-nine 


lead  was  used.  They  conclude  that  no  anomoloua  scattering  is  ob- 
served for  kinetic  energies  y 200  Mev,  and  some  may  exist  for 
lower  energies.  It  is  Interesting  to  note  that  the  curves  for 
0 and  5 cm  lead  show  essentially  the  same  excess  over  their 
theoretical  curves  whether  point  nucleus  scattering  is  Included 
or  not.  It  would  be  desirable  to  have  a direct  measurement  of 
the  momentum  distribution  of  the  incident  particles, 

Leontic  and  Wolfendale  used  a raultlplate  cloud  chamber  at  sea 
level  with  a counter  telescope  that  required  that  detected  par- 
ticles traverse  1 meter  of  lead  above  the  chamber,  six  2 cm  lead 
plates  inside  the  chamber,  and  0,  or  10  cm  lead  below  the 
chamber.  It  was  assumed  that  the  1 meter  of  lead  excluded  all 
but  |x-mesons  from  the  measurements.  The  main  analysis  compared 
the  maximum  scattering  in  any  of  the  center  four  plates  with  the 
r,  m,  s,  angle  for  the  four  plates.  The  analysis  made  the  as- 
sumption that  the  basic  multiple  scattering  law  for  any  given 
particle  should  effectively  be  a gausslan  at  all  angles  if  no 
anomolous  scattering  were  present.  By  an  ingenious  analysis  they 
showed  that  the  results  were  inconsistent  with  this  assumption, 

A further  analysis  assumed  that  an  ’’anomolous  scattering  gaussian” 
would  be  superimposed  on  the  normal  multiple  scattering  gausslan 
in  some  small  fraction  of  the  plate  traversals.  They  then  ob- 
tained a best  matching  of  parameters.  Unfortunately,  the  selection 
criterion  for  considers  1 events  was  biased  in  favor  of  selecting 
cases  where  "anomolous  scattering"  (or  the  non-gaussian  multiple 
scattering  tail)  occured.  Dr,  Wolfendale  Informs  us  that  a 
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preliminary  re-examlnatlon  suggests  that  the  essential  features 
of  the  results  will  not  be  seriously  altered  on  correcting  this 
bias  factor,  but  w«  hare  not  seen  the  details  of  such  a re~examln- 
atlon. 

The  measui^enents  of  Leontlc  and  Wolfendale  have  been  exten- 
ded by  McDlannld^^  ualng  a different  and  Interesting  suialysls, 

20 

Private  ConJinunlcatlon,  We  wish  to  thank  Professor  0,  D, 
Rochester  and  Dr,  A,  W,  Wolfendale  for  correspf'ndlng  with  us 
concerning  the  analysis  of  the  experiments  of  the  Manchester 
group  on  ^-meson  scattering.  The  paper  by  McDlarmid  Is  scheduled 
for  publication  In  Phil,  Mag, 

The  results  seem  to  contradict  those  of  George,  Redding  and  Trent 
in  that  no  anomoloua  fccatterlag  Is  observed  for  low  energies.  At 
higher  energies  the  experimental  results  are  between  the  Mollere 
and  Olbert  multiple  scattering  curves  for  a (partially)  assumed 
distribution  of  Incident  particle  energies.  The  results  are  closer 
to  the  Mollere  than  to  the  Olbert  curves.  We  have  not,  however, 
performed  the  detailed  folding  together  of  our  final  curve  with 
their  derived  momentum  distributions  as  is  necessary  for  a quant- 
itative comparison. 

References  to  earlier  papers  will  be  foimd  in  the  articles 
discussed  ab^ve. 
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FIGURE  CAPTIONS 


Pig.  1,  Curves  illustrating  the  first  method  of  computing  the 

multiple  scattering  distribution  M^(y)  for  cp  = 1 Bev 

in  2 cm  Pb,  including  only  coherent  scattering,  g = 

single  scattering  law*  = multiple  scattering  gaus- 

slan  for  single  scatterings  through  y ^ 0,5;  = 

multiple  scattering  distribution  for  single  scatterings 

through  y > 0,5;  P is  the  assumed  coherent  nuclear 

N 

form  factor.  Multiply  y by  1,74  degrees,  A 

nuclear  size  R = l.OA^'^^  x 10”^^cm  was  used  for  all 
of  the  examples  of  the  figures^ as  discussed  in  Ap- 
pendix A, 

Pig,  2,  Curves  for  the  same  case  as  in  Pig,  1,  The  Moll^re 
and  Olbert  multiple  scattering  distributions  are 
shown  for  comparison  with  M^(y),  The  Olbert  distri- 
bution is  for  a single  scattering  cutoff  angle  of  y « 

1.0,  "The  point  nucleus  euid  extended  nucleus  single 

i 

scattering  distributions  are  g'  emd  g. 

Fig.  3,  Curves  for  cp  = 1 Bev  and  2 cm  Pb  using  for  the 
total  elastic  plus  inelastic  scattering,  g'  = point 
nucleus  single  scattering  law;  g = g*  the  as- 

sumed extended  nucleus  single  scattering  law*  M^ 
is  the  resulting  multiple  scattering  distribution. 
Multiply  y by  1.74  for  Bov  degrees. 

Pig,  4.  Curves  for  cp  = 1 Bev  and  5 om  Pb  using  only  coherent 
scattering,  g^,  g,  and  M^  have  the  same  meaning 


FIGURE  CAPTIONS  (Continued) 


as  in  the  preceding  figures.  Multiply  y by  1.714-  for 
Bev  degrees. 

Fig.  Curves  for  cp  “ 1 Bev  and  5 om  using  the  total  elastic 
plus  inelastic  scattering.  The  symbols  have  the  same 
meaning  as  in  the  preceding  figures. 
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